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Abstract 

Strong-coupling series expansions are calculated for the Hamiltonian version 
of compact lattice electrodynamics in (2+1) dimensions, with 4-component 
fermions. Series are calculated for the ground-state energy per site, the chiral 
condensate, and the masses of 'gluebalF and positronium states. Comparisons 
are made with results obtained by other techniques. 

PACS Indices: ll.Ha .; 12.38Gc. 
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I. INTRODUCTION 



Quantum electrodynamics in 2+1 dimensions (QED 3 ) has generated considerable interest 
over recent years. The model is super-renormalizable, but shares a number of important 
features with quantum chromodynamics (QCD) in 3+1 dimensions: it is believed to be 
confining at large distances (in the quenched approximation, at least), while in the massless 
fermion limit it displays a chiral-like symmetry which is spontaneously broken [|TJ. It is 
thus an ideal laboratory for testing non-perturbative methods of analysis. Versions of the 
model may also be relevant to theories of the new high-T c superconductors f3|. The version 
with two-component massless fermions generates a dynamical mass for the photon through 
a Chern-Simons term Q. This complication can be avoided in the four-component version 
|T]||], where "crural" symmetry is broken in the normal Goldstone fashion, leading to a 
doublet of massless Goldstone bosons analogous to the pion in QCD. For the four- component 
model with Nf flavours of massless fermions, there has been a debate running for some time 
whether chiral symmetry is broken for all values Nf [|1],[§|6| or whether there is a critical 
value N c ~ 3.5 above which no spontaneous symmetry-breaking takes place f/]-||. We shall 
have nothing to say about this question. 

The four-component version has been studied by several different techniques, but we 
remain far from a complete understanding of the model. Euclidean lattice Monte Carlo sim- 



ulations have been performed by several groups [M,^G,11]. A number of authors P,y,|,12,13 



have used Schwinger-Dyson techniques to study the chiral symmetry breaking, and Allen 
and Burden |14]] have also produced estimates of the bound-state meson spectrum at finite 
fermion masses. The Hamiltonian lattice version has been studied by means of strong- 
coupling expansions [15] and a loop expansion technique JR:]. A light-front approach has 



also been discussed ||17yi8|l ; and the non-relativistic limit has been analyzed in some detail 

EM- 

Here we treat the Hamiltonian lattice model by using linked-cluster techniques f2]J to 
generate further strong-coupling series, thus extending the previous results of Burden and 
Hamer [|l5j (hereafter referred to as I). It is well-known that Euclidean Monte Carlo tech- 
niques are difficult and expensive to apply to models with dynamical fermions, and so it 
seems worthwhile to see if other techniques such as strong-coupling expansions can give 
useful information in such cases. A previous analysis of the Schwinger model [2~2"| did indeed 
show that strong-coupling series approximants can converge well into the weak coupling 
region. While not quite as accurate as the exact finite-lattice technique, the series approach 
did give quantitative estimates of the lowest bound-state mass in the continuum limit, at 
about the 5-10% level of accuracy. 

The paper begins with an outline of the lattice formulation of the model in section 
II, followed by a brief summary of the methods of calculation in Section III. Our results 
are presented in Section IV, discussing the ground-state energy, the chiral condensate, the 
"glueball" masses, and the spectrum of the bound-state mesons as a function of the bare 
fermion mass. Our conclusions are summarized in Section V. 



II. FORMALISM 
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A. Continuum Formulation 



The continuum Lagrangian density takes the standard form 

C = ~F^ V + fi-efi-m)^ (1) 

where 

Fp, = d^A u - d u A^ (2) 

and the Lorentz indices //, u = 0, 1 or 2. The electric coupling e in (2+1) dimensions has 
the dimensions of (mass) 1 / 2 . Choosing the timelike axial gauge 



A = (3) 



the Hamiltonian is found to be 



where 



and 



H = J d 2 x{-iip(y + ie + mipip + ^(E 2 + B 2 )} (4) 



E l = F l ° = -A 1 (5a) 



B = d x A 2 - d 2 A l (5b) 

Note that the magnetic field B has only one component in (2+l)D. Here ip is taken as a 
single four-component Dirac spinor [[|, and 70, 71, 72 are 4x4 Dirac matrices, for which we 
shall use the Dirac representation where necessary: 

7o= ( a oJ!)' 7l = l (o-° ffl )' ^(0 -l 2 ) ^ 

In the zero-mass limit, the Hamiltonian (^) possesses a global U(2) "chiral" symmetry 
0, whose Lie algebra is spanned by the matrices 

I, 74= 0) ' 75= *(?" J )' (?) 

and 745 = — ^7475. This symmetry is expected to be spontaneously broken 0, which should 
be manifested by a non-zero value of the chiral condensate (ip^o- 

At large fermion masses, a non-relativistic analysis can be carried out |L8|-[2(J. Cornwall 
and Tiktopoulos and Sen |L9| showed that if the divergences were regulated by giving a mass 
v to the photon, then at one-loop order the renormalized self-mass of the fermion is 

rriR = m + — ml (8) 

47T \ V J 
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while the potential due to one-photon exchange between the electron and positron is: 

V f r ) = _e 2 / d \ = f- L + In —) + 0(v 2 r 2 ) (9) 

V ; J (2tt) 2 k 2 + u 2 2tc\' 2 J K J K J 

where 7 is Euler's constant. Both quantities show logarithmic divergences, but these diver- 
gences cancel in the Schrodinger equation for the positronium bound states 

1 - 1 - e 2 

EVOr) = [ V 2 , + 2(m R - m) + V(r)Wr) = [ Vl + —(7 + \nmr)U(r) (10) 

m m27T 

Numerical solutions of this equation |T8^pO| give the 'binding energies' of the lowest positro- 
nium states as 

ES = 1^.7968- ilng)) (Ha) 

£?= £ (2.9323- \ ln(gj (lib) 

for angular momentum Z = 0, and 

B -£l( 2 ,5 6 6-i,„g!)) (12, 

for 1 = 1. At leading order the binding energies are independent of "spin", so that each of 
these energy levels should be four-fold degenerate in the four-component fermion model. 



B. Lattice Formulation 



A 'staggered' Hamiltonian lattice formulation of this model has been discussed in refer- 
ence I. The four components of the continuum fermion field fit naturally onto the four sites 
of a 2 x 2 unit cell on the 2-dimensional spatial lattice, leading to a lattice Hamiltonian as 
follows: 

H = fw (13) 

where 

W = W + yW 1 + y 2 W 2 (14) 

andQ 



The fermion mass term given in reference I had the wrong sign. As it turns out, that did not 
affect the results for the quantities they calculated. 
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Wo= We + Wii = STE? + ^£(-ir+ r2+ Y(r-Mf) 

l r 

Wi= ^Vi(f)[x ] ~(f)Ui(f)x(r '+ i) +h.c] 

r,i 



(15a) 
(15b) 
(15c) 



Here f — (ri,r 2 ) labels the sites, / the links, p the plaquettes and i = 1,2 the directions on 
a square two-dimensional spatial lattice with spacing a. The dimensionless coupling g and 
mass parameter fi are defined in terms of their continuum counterparts e and m by 



g 



e a and /i 



2am 2m 
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(16) 



while y = 1/g 2 , and ?7i(r) = (— l) r2+1 , 772 (r) = 1- The term VF e is the electric field term, 
Wfj, is the fermion mass term, W± is the fermion kinetic energy, and W 2 is the magnetic field 
energy, involving the usual plaquette operator U p . 

The correspondences between the lattice fields and their continuum counterparts are for 
gauge fields 



a 



A 



E\x) 
A\x) 



where the link operator 



Ui = exp[ieaAi(r) 
while for the fermion field components 



-i 1 

1 -i 
-i 1 
1 -i 



(17a) 
(17b) 

(18) 

























^3 




.6. 




.V"4. 



(19) 



where 23 



£(0 



•ri+r 2 



x(0 



(20) 



and the components 1, • • • , 4 are assigned to sites of the 2x2 unit cell as shown in Figure 1. 
The commutation relations between the lattice fields are 



[E l ,U l i\ = U l 5 w 
[E ll Ul\ = -U}5 lv 

{xHr),x(7)} = kr- 



(21a) 
(21b) 
(21c) 
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[E h x(r)] = [E h ^(r)] = P h X(f)] = [U h x f (f)] = (21d) 

With these correspondences, it can be shown |15| that the lattice Hamiltonian (|T3|) reduces 
to the continuum Hamiltonian @) in the naive continuum limit a — > 0. 

The introduction of the lattice breaks the U(2) "chiral" symmetry down to a discrete 
symmetry generated by shifts of one lattice spacing jl5| . A unit shift in either the x or y 
direction leaves the kinetic term in the Hamiltonian (p~3[ ) invariant, but alters the sign of the 
mass term. The corresponding continuum field transformations, from (19), are 

V> -> e i(7r/2)74 ^ (22a) 

or 

jj, _> e^/ 2 ^ (22b) 

respectively. 



C. The Strong-Coupling Limit 

The Hamiltonian ( |T3"D acts an a Fock space spanned by the usual strong-coupling basis 
24]. With each link is associated an integer electric flux n\ such that Ei\ni) = ni\ni). The 
operators Ui and U} increase and decrease the flux on link / by one unit, respectively. Each 
site of the lattice can be in one of two fermionic states |+) or |— ) obeying 

X j H = \+}, X f |+) = (23a) 
x|-) = 0, xl+> = |-> (23b) 

Consider first the massless theory, fi = 0. In the strong-coupling limit, the variable y = 
and the Hamiltonian W reduces to W e . The ground-state is then highly degenerate, having 
flux m = on each link, but with the fermionic state entirely arbitrary. This degeneracy 
is broken at the next order by the kinetic term W±, leaving only two degenerate states \A) 
and \B) whose fermionic content is 

i A \ fl+)> on odd sites; ^ 
1 1— ), on even sites. 

and 

N-), on odd sites; 
L |+), on even sites. 

The chiral shifts of equation ([22|) map these two states into each other. When the mass term 
Wfj, is included, chiral symmetry is explicitly broken and state \B) is favoured energetically. 
We thus take \B) as the unperturbed strong-coupling ground state for both the massive and 
massless cases, and interpret this as the state with no fermion excitations present. 

An excitation on an odd or even site creates a positively or negatively charged fermion 
respectively, i.e. a positron or electron. The first-order perturbation W\ creates or destroys 
an electron-positron pair on neighbouring sites, joined by a link of flux. The second-order 
perturbation term Wi creates or destroys a plaquette of flux. Gauge invariance ensures that 
for any state obtained from the unperturbed vacuum by application of the operators W\ and 
W2, the net flux from any site is equal to the charge of the fermion at that site, i.e., Gauss' 
law is obeyed. 
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D. Positronium States 



This theory is expected to display confinement [IS], and the only fermionic states with 
finite energy are expected to be electrically neutral "positronium" bound states. In the 
strong-coupling limit, the lowest energy positronium states consist of an electron-positron 
pair on neighbouring sites, connected by a link of unit flux. There are eight translationally- 
invariant states of this type, corresponding to the eight links in the unit cell, and we need to 
identify the linear combinations of these states which correspond to eigenstates of the lattice 
symmetry operators. The corresponding procedure for meson states in four- dimensional 
Euclidean lattice QCD has been discussed by Golterman [^Sj . 



The symmetry group of the lattice Hamiltonian ( |13l) is composed of the following ele- 
ments: 



1. Even translations 
xif) -> xif + 2i), U^r) -> U j (r + 2i) (26) 
This corresponds to spatial translational invariance in the continuum model. 



2. Odd translations 
X(r)-x(r-+i), Uiifi^Uiif+l) (27a) 

or 

X (f) - (-l) ri x(r + 2), Ui{r) - U i {r + 2) (27b) 

These are the discrete lattice versions of "chiral" symmetry corresponding to equations (p2|) . 
The massless Hamiltonian is symmetric under these transformations, but not the massive 
one. 



3. Diagonal shift D 



A combination of a shift by one site in the x direction and one site in the y direction 
gives a diagonal shift 

X (r) - (-l) ri x(r + 1 + 2), U t (r) ^ U l (f+ 1 + 2) (28) 

This corresponds to a discrete 745 rotation in the continuum fields 

i> -> ^745^ (29) 
This remains a symmetry of the massive Hamiltonian also. 
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4- Square lattice rotations, R 

Let R denote a lattice rotation by 7r/2 about a perpendicular shown by Fig. 2: 

X(f) -> ^(r')x(r') (30a) 

C/ 2 (r) -> f/i(f') (30b) 

-» f4(r'-2) (30c) 

where 

ri = r 2 , r 2 = -n (31) 

and 

i2(r x ,r 2 ) = -[(-1)^ + (-lp + ("l) ri+r2 " 1] (32) 

Repeated rotations generate the rotational symmetry group of a square, with 4 elements. It 
corresponds to rotation in both space and "spin" in the continuum model. 

5. "Axial parity" inversion, A 



The "axial parity" inversion is discussed by Burden and Allen ]T2,H|. In the continuum, 
it corresponds to the operations: 

ip(x) -> VV) = i>{x) -> ip'(x') = i){x)A~ 1 (33) 

and the vector field transforms as: 

A(x) -> A'(x') = -i*(x) (34) 

where ). A suitable representation for the fermion operator A is the 

matrix Z7 . 

On the lattice, this is simply 

Xif) -> (35a) 
C/i(f) -> U}(-f-i) (35b) 

(35c) 

which is equivalent to i? 2 , a rotation by II in 2+1 dimensions. 

6. Reflection, II 

A reflection in the y axis corresponds to 

Xif) -> x(0 (36a) 

^(r) -> C/}(F'-i) (36b) 

t/ 2 (r) -> C/ 2 (f') (36c) 

where 

r[ = — ri, r 2 = — r 2 (37) 



S 



7. Charge conjugation, C 



The charge conjugation operation is also discussed by Burden and Allen [T^JT^]. In the 
continuum, it corresponds to the operators: 

ip(x) -> tp'(x) = C^ T (x), tp(x) -> ij'(x) = -ip T (x)C^ (38) 

and 

A»{x) -> A^'(x) = -A"(ar) (39) 

The fermion operator can be represented, up to an arbitrary phase, by the matrix 7 2 . The 
translation to lattice variables is slightly involved, and is best handled in terms of the fields 
£i defined in equations (|T9|), (|20|) . We shall not go into further details here. 

The translationally invariant positronium eigenstates can be classified in terms of their 
eigenvalues under these symmetry operations. The group of square rotations is C4, the cyclic 
group of 4 elements. It has 4 irreducible representations, each with dimension 1. The allowed 
eigenvalues of R simply consist of the 4th. roots of unity (i.e. powers of e = e~~ 27 ™/ 4 ). The 
shift D, reflection II and axial parity A each generate 2-element groups, with eigenvalues 
±1 for the positronium states with translational symmetry. 

In the strong-coupling limit, the "link excitations" on the unit cell corresponding to 
the lowest energy positronium states consist of the eight states listed in Table 1, numbered 
according to Figure 3. These states transform into each other under the action of the 
symmetry operators. By taking linear combinations of these states, one can form eigenstates 
\ipi), • • ■ , \ips) of the lattice symmetry operators, as listed in Table 2. Note that the states 
with rotation eigenvalue R = e = — i or R = e 3 = — i cannot simultaneously be eigenstates 
of II, because a reflection converts R = i to R = —i, and vice versa. We have chosen 
to list states to \ip$) which correspond to eigenstates of II, and are thus symmetric 
or antisymmetric combinations of the states with R = ±i. Similarly, we have chosen to 
list states \ip2) and ^4) which are eigenstates of C, rather than R. All eight states are 
degenerate in energy in the strong-coupling limit y — > 0. The fact that the Hamiltonian 
is symmetric under both rotations and reflections implies that the pair ^5) and l^) will 
remain degenerate at all couplings, and likewise the pair ^7) and \if>s). The combination 
of rotation and charge conjugation symmetry implies that the pair an d IV^) will also 
remain degenerate. 

In the 'naive continuum limit' a — > 0, when U\ — > 1 + ieaAi, one finds that the quartet 
of states I Vs) to iV's) transcribe to simple combinations of quark and antiquark fields on the 
lattice, and correspond to 'vector' states in the language of Burden and Allen p!2]JT4| , with 



J AC = 1 , where J is the 'total angular momentum'. The states \ip\) to ^4), on the other 
hand, contain an admixture of gauge fields in the naive continuum limit, and have no direct 
counterparts in the catalogue of quark-antiquark states discussed by Burden and Allen. The 
state is a scalar state, having the same quantum numbers as the vacuum. 



E. Weak-coupling expansion 

Some useful information on the ground-state properties, at least, can be gained by per- 
forming a "weak-coupling" expansion for the lattice system as t/ — > 00. For the present 



9 



model, this exercise was carried out in reference I. 

The ground-state energy per site has an asymptotic expansion in the weak-coupling limit: 

u ~ -2y 2 + 1.9162?/ / dq x / cWcos 2 q\ + cos 2 q 2 + -r^V + Oil) as y —> oo 

n z Jo Jo Ay z 

(40) 

where the integral arises from a sum over the fermionic degrees of freedom. At very large y 
one finds 

coo 2y 2 + 0.9581?/ + 0(1) as y -> oo (41) 

but for finite y and /i it is more useful to evaluate the integral as it stands. 

The ground-state expectation value of the chiral condensate was also calculated [T5| as 

(^) lattice y- dq x A/ 2 [cos 2 Ql + cos 2 q 2 + ^Y 1 ' 2 + 0(y- 2 ) as y - oo 

-K l y Jo Jo 1y l 

(42) 

where 

(^) lattice = (H^i-iP^xHmmo) (43) 

r 

= continuum ^ 44 ) 

in terms of the continuum chiral condensate. This quantity is given by the Feynman Hell- 
mann theorem as 

= ^ (45) 
where u is the ground-state eigenvalue of W. 



F. Non-Relativistic Limit m/e 2 — ► cxd 

When the fermion mass m becomes very large, it should be possible to study the model 
in a "quenched" approximation, where fermion loop diagrams are suppressed. In a staggered 
lattice formulation such as the present one, the suppression of all fermion loops will lead 
to the "static" limit, in which any fermion excitation is fixed at its initial lattice site, and 
apart from the mass term the remaining lattice Hamiltonian is simply that of the pure gauge 
theory. 

To allow the fermions to move or migrate on the lattice, one has to go to the next order in 
powers of l//i, the mass parameter, and allow second-order diagrams involving the excitation 
of an (e + e~) pair on neighbouring sites, as shown in Figure 4. Figure 4a) represents a loop 
diagram in the vacuum sector, and 4b) a similar diagram in the one-fermion sector; Figure 
4c) illustrates a "hopping" diagram, in which an existing fermion hops two lattice sites. This 
allows fermion migration to take place. 



10 



It is useful to define new fermion variables representing the excitations on the strong- 
coupling ground-state, namely: 

m = { X ^^ ( r i + r 2)even; ( } 

\x(r), (n + r 2 ) odd. v ; 



so that 



0(r)|O) = O (47) 



where |0) = \B) is the strong-coupling ground state. The Hamiltonian ([TJ]) and link excita- 
tions (Table 1) are easily translated in terms of the new variables. 

In the quenched approximation outlined above, the effective Hamiltonian in the vacuum 
sector is 

AT ,,2 

= -N y - + W e + y 2 W 2 (48) 

(where N is the number of sites on the lattice) which is simply the pure gauge field Hamil- 
tonian plus a constant. The first constant term is the negative energy of the "Dirac sea" on 
the lattice, while the second constant term is the contribution of the second-order diagram, 
Figure 4a). More interesting is the effective Hamiltonian in the (e + e~) sector, which takes 
the form: 

Weff +e ~ = - (iY ~ 4 V - (N - 4) V - + W e + y 2 W 2 + y 2 W[ (49) 



2 r y ' /i 



where 



W[ = — (-lf m {^(f- i)U\(r- l)[Ul(r)(j)(f+ 1) + (-l) r2+1 ([/ 2 (f - 2)0(r - 2) 

f odd 

+U\(r)<P(r + 2))] + <j,\r - 2)U ] 2 (r- 2)[U 2 {r)<f){f + 2) + (-l) r2+1 [/ 1 t (r)</)(f + 1)] 
+(-l) ra+ V t (r + 2)?7 2 (r)?7i(r)0(f + 1) + h.c.} 

+ ^ E (-^^{^(r-^^i^-^^^f+i) 

' r even 

+(-l) r2+1 (C/ 2 f (r - 2)0(f - 2) + U 2 {f)<f>{f + 2)] 

+</> t (f- 2)C/ 2 (f- 2)[f/ 2 f (f)0(f + 2) + (-l) r2+1 [/ x (r)</>(r + i)] 

+ (-l) r2+1 t (r + 2)^(70^1(0^+ 1) + hx -} ( 50 ) 

This clumsy expression is merely a constant, plus the pure gauge Hamiltonian, plus 
"hopping" terms for the fermions in the six different paths allowed for a double hop on the 
staggered lattice. The associated phase factor 

f^Hi?) = { + 1 ' if site ^occupied; > 5 _ 
~ \ — 1, if site r unoccupied. 

accounts for the change of sign if the hopping fermion "passes through" an occupied site. 
These hopping terms correspond to the kinetic energy in the non-relativistic continuum 
Hamiltonian. 
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The effective Hamiltonian is more complicated in form than the original lattice Hamilto- 
nian, but does not allow any further fermion excitations, and would therefore be quicker and 
easier to implement in numerical calculations. We have not attempted any such calculations 
as yet. 



III. METHOD 

To calculate the strong-coupling series for the model, we used Nickel's cluster expansion 
method. The techniques necessary were reviewed in He et al [21|, and will not be repeated 



here. In these calculations, the Wq in Eq. ([14]) is taken as the unperturbed Hamiltonian, 
diagonal in the basis of eigenvectors of E h while the W\ and W 2 in Eq. flUf ) then act as 
perturbations. 

To generate the series for the ground state energy, we need to generate a list of connected 
plaquette configurations, together with their lattice constants and embedding constants. 
Since the first-order perturbation W\ and second-order perturbation W2 involve links and 
plaquettes, respectively, a cluster a will contribute terms 0(y a ), where a is given by 

a > 2n p + m (52) 

where n p is the number of plaquettes in a, and n\ the number of links not contained in 
plaquettes. Up to the order y 22 considered in the current paper, there is only one graph 
(Figure 5) which does not obey the above relation, it actually contributes to order y 22 (due 
to the combination of W2 on each plaquette and W\ in outer links) instead of y 24 according 
to Eq . (|52"D . There are a total of 5494 graphs which contribute up to order y 22 for the 
ground-state properties. 

The calculation of glueball masses involves a list of clusters, both connected and discon- 
nected, with at least one plaquette in each graph. There are 457 graphs which contribute 
to order y 10 . 

The calculation of meson masses generally involves a list of both connected and discon- 



nected clusters [21]. The eight different links in the unit cell (shown in Figure 3) are not 
equivalent in those calculations, which means we cannot identify clusters which are topo- 
logically equivalent, or even use rotation or reflection symmetry. Thus the separate clusters 
proliferate enormously in this case: there are 164142 clusters contributing up to order y 10 . 
For the scalar meson mass mi, the eight different bond types are equivalent, so we only need 
actually one bond type, and there are only 569 clusters which contribute up to order y 12 . 



IV. RESULTS 

A. Ground-state energy 

Using the linked-cluster expansion method, series have been calculated for the ground- 
state energy per site up to order y 22 . The first few terms are: 

uj /N= — - - — + 



2 l + 2/i 2 (l + 2//) E 
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y 6 (-4742 - 5084/t - 1640/i 2 - 368/i 3 - 64/i 4 ) 
(1 + 2^)5(3 + 2/i) (7 + 2/i) 



(53) 



where ujo is the ground-state eigenvalue of W. The coefficients are listed for various fixed 
values of the dimensionless mass parameter fi = 2m/ e 2 in Table III. These coefficients agree 
with those of reference I up to 0(y 16 ). 

Extrapolating these series into the weak-coupling region using integrated differential 
approximants [Q, one obtains results as shown in Figure 6. For the large-mass case, /i = 10, 
it can be seen that the strong-coupling approximants match on to the weak-coupling form 
(fi(i|) very nicely at around 1/y — 0.5. For the lower masses, the strong-coupling approximants 
do not converge well enough to establish a precise matching, but they are clearly quite 
consistent with the asymptotic form (]40f) . It is noteworthy that the //.-independent form 
( PD is only attained for very weak couplings /i 2 /y 2 < 0.4. 

The successful matching between the strong-coupling approximants and the weak- 
coupling asymptotic form gives some confidence that the series coefficients have been calcu- 
lated correctly, and that the approximants converge well enough to provide useful informa- 
tion about the weak-coupling (continuum) behaviour. 



B. Chiral condensate 



A more interesting quantity is the chiral condensate. The first few terms in the series 
are: 

X ^ 7 2 1 + 2// (1 + 2/i) 4 

% 6 (122987 + 245156/t + 181668/i 2 + 64656/x 3 + 13392/i 4 + 2048/i 5 + 192/i 6 ) 
+ (l + 2//) 6 (3 + 2/x) 2 (7 + 2/i) 2 + ' 

Further coefficients at fixed values of \x are listed in Table III. Figure 7 shows the extrapola- 
tion of these series into the weak-coupling region, as compared with the weak-coupling form 
(fj2f). Once again, the large mass results (/t = 10) match very well to the weak-coupling form, 
while the lower-mass ones have not yet attained it before convergence of the strong-coupling 
approximants is lost. 

The most interesting case is /t = 0, the zero mass limit. A graph of y 2 (ipif)) la,ttlce against 
1/y is shown in Figure 8. It can be seen that integrated differential approximants to the 
series fail to converge below about 1/y — 1.5, although Pade approximants behave in a more 
consistent fashion. If taken at face value, the Pade approximants would indicate a very large 
value of the chiral condensate in the continuum limit l/j/ — ^ 0, 

e~ 4 (^) physical = -0.284(10) , (55) 

Also shown in Figure 8, however, are some Monte Carlo estimates of the chiral condensate 
for the Euclidean version of this modelR by Burkitt and Irving pT|. Their results are roughly 



2 In making this comparison, we assumed the Euclidean coupling gE and Hamiltonian coupling 
gn are equal, lacking information on a more precise connection. 
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compatible with ours at about 1/y — 0.6, but show a dramatic decrease in magnitude beyond 
that point. Since this occurs well below the region of convergence of the series approximants, 
the series provide little evidence either to confirm or deny this phenomenon. 

Another Monte Carlo calculation of the chiral condensate in the non-compact version 
of the model has been carried out by Dagotto, Kogut and Kocic |7[|. They found a quite 
different behaviour, in which f3 2 (ipijj) lattlce (where (3 = y = 1/g 2 ) plunges rapidly toward zero 
at a low P value around (3 ~ 0.4, and then levels out to a plateau at a very small value, 
around j3 2 (ipip) ~ 0.001. Although this value is only tracked to (3 ~ 1, they take this as 
evidence of chiral symmetry breaking in the continuum limit. Burden and Roberts have 
also obtained a Schwinger-Dyson estimate of the chiral order parameter, 

e" 4 (^) ~ -0.003 (56) 

which is at least a similar order of magnitude to that of Dagotto et al. [|7|. It differs by a 
full two orders of magnitude from the naive extrapolation above, equation (|55|) . 

It would be interesting to see a more extensive Monte Carlo simulation of the compact 
lattice model, to check whether the decrease seen by Burkitt and Irving JXTJ] is real, and 



whether the chiral condensate subsequently levels out at a small plateau value as in the 
noncompact model. There seems to be very little prospect that further series calculations 
could shed light on these questions. 



C. 'Glueball' masses 

Strong-coupling series for the 'photonball' masses, rriA and ms, corresponding in the 
strong-coupling limit to single plaquette excitations which are antisymmetric and symmetric 
under reflections, respectively, have been calculated to order y 10 . The leading terms for both 
series are: 

m SA = 4 + 16y 2 /((l + 2/i)(l - 2/i)(3 + 2/i)) + 0{y 4 ) (57) 

with the difference between the two series only emerging at order y 4 . Further coefficients at 
fixed values of \i are listed in Table IV. The coefficients up to 0(y s ) were previous calculated 
in reference I. 

Approximants to these strong-coupling series are graphed in Figures 9 and 10, and com- 
pared with the results for the pure gauge theory, which we have calculated previously [^7] to 



order y 32 . It can be seen that for large /i values the series behave very similarly to the pure 
gauge case, and appear to be decreasing exponentially towards zero as y increases. Our re- 
sults are consistent with those of Burkitt and Irving [|TTJ, who found a systematic downward 



shift in rriA for moderate quark masses, as seen for the case /x = 2 in Figs. 9 and 10. At 
small \i values, however, the behaviour appears to change somewhat. In the region /j ;$ 1, 
there is substantial overlap and level crossing between the glueball and meson states in the 
strong-coupling region, and the picture becomes more confused. The series approximants 
show an initial rise for the glueball masses in the strong-coupling region, followed by an 
apparent turnover, at least in the /i = 0.5 case, but the convergence is not sufficient to track 
the behaviour reliably into the weak-coupling region. 
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D. Positronium masses 



Strong-coupling series for the meson states discussed in Section II have been calculated 
to 0(y 10 ) (or to 0(y 12 ) for mi). The first three terms for arbitrary // are 

m 1 = l + 2n + Uy 2 /(l + 2/i) - y 4 (535 + 186// + 60/i 2 + 8// 3 )/[3(l + 2/i) 3 ] + 0(y 6 ) (58a) 

m 2 = m4 = 1 + 2/i + 10y 2 /(l + 2/i) - y 4 (283 + 42/i + 12/i 2 + 8/i 3 )/[3(l + 2/i) 3 ] + 0(y % ) (58b) 

m 3 = 1 + 2/i + 6?/7(i + 2/i) + y 4 (-175 + 6/i + 36/i 2 - 8/i 3 )/[3(l + 2/i) 3 ] + 0(y 6 ) (58c) 

m 5 = m 6 = 1 + 2/i + 6i/ 2 /(l + 2//) - i/ 4 (199 + 66// + 12// 2 + 8// 3 )/[3(l + 2//) 3 ] + 0(i/ 6 ) (58d) 

m 7 = m 8 = l + 2/i + 6i/ 2 /(l + 2//) - y 4 (199 + 66// + 12// 2 + 8// 3 )/[3(l + 2//) 3 ] + 0(i/ 6 ) (58e) 

Further terms at selected values of /i are given in Table V. 

Figure 11 graphs the various meson masses as functions of 1/y, at a fixed, large mass 
parameter /i = 10. It can be seen that the series approximants for these masses do not 
converge below about 1/y — 1.5. For the lowest vector meson states there is a suggestion 
that the mass reaches a peak at 1/y ~ 1.5, and then turns downward. A crude linear 
extrapolation has been made to estimate the continuum limit, but the uncertainly in the 
estimate is large. 

The resulting values for the vector mass are graphed as a function of // in Figure 12, 



along with the Schwinger-Dyson estimates of Allen and Burden [14], and the non-relativistic 



prediction Hi~8|-f20|j. It can be seen that the lattice estimates of the positronium binding 



energy have large errors, especially at smaller /i, and lie about three times higher than 
the Schwinger-Dyson values. Neither the Schwinger-Dyson nor the lattice estimates show 
any definite evidence of the logarithmic increase at large /i predicted by the non-relativistic 
theory. 

Figure 13 shows the masses as functions of 1/y for the massless case /i = 0. Once again, 
convergence is lost at rather small y values, around 1/y — 2.5, and it is hardly possible to 
make useful estimates of the continuum limit. There is no sign of any of the masses dropping 
towards zero, or acting like a Goldstone boson. This is because the expected Goldstone 



boson states are the "axiscalar" and "axipseudoscalar" states ||14j| , which are not among the 
single-link excitations in the strong-coupling limit which we have treated here (see Section 
IID). The Goldstone bosons probably correspond to L-shaped double-link excitations, which 
will transform into each other under a single plaquette excitation. It would be interesting 
to study their behaviour, but we have not yet attempted such a study, owing to technical 
complications. 



V. SUMMARY AND DISCUSSION 

New strong-coupling series have been presented for the ground-state energy and chiral 
condensate, along with the 'glueball' and positronium masses, in 4-component Hamiltonian 
lattice QED2+1 with full dynamical fermions. This represents the first attempt at a series 
calculation for the positronium states in this model. 

Two major features are evident from these results. Firstly, there is a very clear sepa- 
ration of scales between the 'glueballs' and the positronium states as the continuum limit 
is approached. At large fermion mass /i, the positronium energies remain finite in the 
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continuum limit, while the 'glueball' masses scale exponentially towards zero, presumably 
corresponding in the limit to massless photon states. The same thing appears to happen at 
smaller /i values, although our evidence for the exponential decrease of the glueball mass is 
rather slim. This separation of scales is consistent with earlier discussions of the pure gauge 



model |]28fl : in the nonrelativistic or static fermion limit, if one sticks to naive 'engineering' 
dimension scales one will end up with a theory of free, massless photons as above, whereas 
if one renormalizes the scale as discussed by Polyakov [^] or Gopfert and Mack [|3U] one will 
obtain a theory of free, massive bosons. 

The second feature is that the addition of dynamical fermions does not greatly affect the 
glueball masses at large fi: at \i = 2, for instance, the only effect was a reduction of the 
glueball mass at fixed y of order a percent or two. For % 0.5 the effect is more pronounced, 
however. 

In general, the results of these calculations were somewhat disappointing. The bulk 
ground-state energy per site converged well enough into the weak-coupling region to display 
a convincing match with analytic weak-coupling expansions jL5|, and the 'glueball' masses 
converged well enough to justify the statements given above. Other quantities, however, 
were not mapped out with such success. Series approximants to the chiral condensate at 
H = only converged down to l/y ~ 1.5, and were unable to confirm or deny the rapid 
plunge in magnitude seen by Burkitt and Irving [] 1 around l/y — 0.5. It will require more 
detailed Monte Carlo studies to confirm that this plunge really occurs, and to determine 
whether the chiral condensate then levels out at a small but finite value. 

Series approximants to the positronium masses likewise only converged down to l/y — 
1.5, even at large mass fi, making any extrapolations to the continuum limit very uncertain. 
Our estimates of the continuum 'binding energy' are thus of little more than qualitative 
accuracy. They lie about three times higher than the Schwinger-Dyson estimates of Allen 
and Burden [[14[], which are not expected to be very accurate at large mass /i in any case. 
Neither set of results shows any definite sign of the logarithmic increase in energy with /x 
predicted by the non-relativistic analyses |18|-f20"|| . Once again, a detailed understanding 



of the positronium spectrum must await a more accurate study, by Monte Carlo or other 
methods. 

The question whether this model develops Goldstone bosons in the massless limit fi = 
was not explored. These 'axiscalar' and 'axipseudoscalar' states are not among the 
single-link excitations for which we have calculated strong-coupling expansions. They are 
most likely to be found among the 'double-link' excitations. It is technically more difficult 
to calculate strong-coupling series for the double-link excitations; but on the other hand, the 
series may well converge more quickly for these low-lying excitations. This might provide 
an interesting subject for further study. 

In general, however, the prospects for further series calculations look dim. In contrast 
to the Schwinger model case |[22|| , the strong-coupling series approximants do not converge 
far enough into the weak-coupling region to allow accurate extrapolations to the continuum 
limit. Thus it seems there is still no real alternative to Monte Carlo methods, whatever 
their limitations, for the detailed investigation of lattice gauge theories in three and four 
dimensions. The strong-coupling series provide an accurate 'platform' of results in the 
strong-coupling region, which can provide useful calibration points for other methods; but 
their extrapolation into the weak-coupling region remains rather uncertain. 
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FIGURES 




(n,r 2 ) = (0,0) 

FIG. 1. Assignment of spinor components to sites of the 2x2 unit cell. 




(n,r 2 ) = (0,0) 

FIG. 2. A square lattice rotation by ir/2. Point 1 is assumed fixed. 
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FIG. 3. Numbering of sites and "link" states on the unit cell. 
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a) b) c) 

FIG. 4. Lattice perturbation theory diagrams at second order involving excitation of an (e + e~) 
pair. The associated gauge field excitations are not shown. Horizontal lines represent the action 
of 'link' excitation operators from W\\ vertical lines represent the resulting fermion excitations. 



FIG. 5. The last graphs contributing at order y 22 . 

FIG. 6. Graph of the ground-state energy per site, y~ 2 u>o/N versus 1/y, for various fixed 
values of the mass parameter [i = 0,0.5,2,10. The curves at large 1/y are integrated differential 
approximants to the strong-coupling series; while those at small 1 jy correspond to the asymptotic 



weak-coupling form (40) 



FIG. 7. Graph of the chiral condensate, ^-^ lattlce versus [i/y, for various finite values of 
the mass parameter = 0.5,2,10. The curves at large 1/y are various integrated differential 
approximants and Pade approximants to strong-coupling series, while the solid line at small ji/y 
is the weak-coupling asymptotic form (fi"2|), 

FIG. 8. Graph of the chiral condensate, y 2 {"i[}ij))^ ttlce versus 1/y, for the massless case /i = 0. 
The curves shown are integrated differential approximants (solid lines) and [n/(n + 1)] Pade ap- 
proximants (broken lines) to the strong-coupling series. The last two Pade approximants are almost 
indistinguishable, and the successive intercepts at 1/y = are 0.232,0.274,0.283 and 0.284. 

FIG. 9. The mass rriA as a function of y, for [i = 0,0.5,2,10, Various integrated differential 
approximants and Pade approximants to strong-coupling series are shown, together with corre- 
sponding results for the 'pure gauge' case (// = do). 
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FIG. 10. The mass mj as a function of y, as in Figure 9. 

FIG. 11. Strong-coupling series approximants to m//x as functions of 1/y for the meson states 
mi, • • • , mg, at mass parameter fi = 10. 

FIG. 12. Graph of m//x versus /U for the vector state m-j. 

FIG. 13. Series approximants to mi, • • • , mg versus 1/y at /x = 0. 
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TABLES 



TABLE I. Link excitations on the unit cell corresponding to positronium states in the 
strong-coupling limit, with sites labelled according to Figure 3. 



Link state Operator equivalent 



11) 


X t (2)C/ 1 T (l) X (l) 


|2> 


X t (2)f/ 1 (2) X (5) 
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X t(4)C/ 2 t (l) x (l) 
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X t(4)C/ 2 (4) X (8) 


15) 


xt(4)EM4)x(3) 
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X t (6)t/ 1 t (3)x(3) 
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X t (2)[/ 2 (2)x(3) 


18) 


X t (7)t/ 2 t (3)x(3) 



TABLE II. Linear combinations of the link states forming eigenstates of the lattice symmetry 



group = a>i\i}), and their symmetry eigenvalues. 
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TABLE III. Series coefficients of y 2n in strong-coupling expansions of the ground-state energy 



loq and chiral condensate (tpi/i) 
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-1.327976394916X10 7 


-1.291627245890X10 4 


-2.768670289809 




-9.300630428741x10" 


-7 


10 


8.584509655365 xlO 13 


1.206675388588 x10 s 


5.492904223987X10 4 


4.708649820309 




1.110110045268x10" 


-7 


11 


-3.000444314331 xlO 15 


-1.107246112597X10 9 


-2.359089145008 xlO 5 


-8.090428010991 




7.857720150141x10- 


-7 



TABLE IV. Series coefficients of y 2n in strong-coupling expansions of the symmetric and anti- 
symmetric glueball mass gaps 777,5 an d 777^4 of the dimensionless Hamiltonian W. 



n = 



: 0.5 



1 



= 2 



10 



symmetric glueball mass gaps m$ 



4.000000000000 

1 5.333333333333 

2 -4.189814814815X10 1 

3 5.935307344636 xlO 2 

4 -3.229319000985 xlO 3 

5 -6.430232339023 x10 s 



4.000000000000 

1.500000000000 
-3.135416666667 
-1.170910493827xl0~ 2 

8.663828760713X10 1 
-1.617503042067X10 3 



4.000000000000 
-1.066666666667 

4.769703703704 
-2.394715197442X10 1 

1.501264118830X10 2 
-8.784848865979 xlO 2 



4.000000000000 
-1.523809523810X10- 1 
-2.352469495735X10" 1 
-4.890552931893X10- 1 

2.073466005647 
-5.575125327898 



4.000000000000 
-1.743489157677xl0~ 3 
-4.165631870342X10" 1 
-3.937708897154 xlO" 4 
4.945486284578X10- 1 
1.071338323552xlQ- 3 



antisymmetric glueball mass gaps m,A 



4.000000000000 

1 5.333333333333 

2 -4.306481481481 xlO 1 

3 6.055105324434 xlO 2 

4 9.053923337359 xlO 3 

5 -1.377596830177X10 6 



4.000000000000 

1.500000000000 
-4.302083333333 

1.838348765432xl0~ 2 

5.354226523059 xlO 1 
-5.970881550866X10 2 



4.000000000000 
-1.066666666667 

3.603037037037 
-2.776504371231 xlO 1 

1.574287286646 xlO 2 
-9.103588055806X10 2 



4.000000000000 
-1.523809523810X10- 1 
-1.401913616240 
-3.473208512644X10- 1 

1.695152914866 
-7.260136671150X10- 1 



4.000000000000 
-1.743489157677X10- 3 
-1.583229853701 
-1.358230163160X10- 3 

1.162208030611 

2.449261377761 xlO" 3 
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TABLE V. Series coefficients of y 2n in strong-coupling expansions of the meson mass gaps rrii 
(i = 1, 2, • • • , 8) of the dimensionless Hamiltonian W. 



mi 



n%2 = rri4 



m 3 



m 5 = m 6 



mj = mg 



1.000000000000 

1 1.400000000000 xlO 1 

2 -1.783333333333X10 2 

3 3.625151675485 xlO 3 

4 -9.101099254509X10 4 

5 2.534993350264 x10 s 

6 -7.521807853058 xlO 7 



1.000000000000 

1.000000000000 xlO 1 
-9.433333333333 xlO 1 

1.757293650794X10 3 
-4. 164808339711 xlO 4 

1.113036398112X10 6 



(i = 
1.000000000000 
6.000000000000 

-5.833333333333 xlO 1 
1.200151675485 xlO 3 

-2.976672075609 xlO 4 
8.177458288471 xlO 5 



1.000000000000 

6.000000000000 
-6.633333333333X10 1 

1.341849206349X10 3 
-3.324256060799 xlO 4 

9.143840136437X10 5 



1.000000000000 

6.000000000000 
-6.633333333333 xlO 1 

1.253207231041 xlO 3 
-3.008509651089 xlO 4 

8.115398963424X10 5 



2.000000000000 

1 7.000000000000 

2 -2.683333333333 xlO 1 

3 1.365322916667xl0 2 

4 -8.920838363922X10 2 

5 6.436060668280 xlO 3 

6 -4.984333949987X10 4 



2.000000000000 

5.000000000000 
-1.283333333333X10 1 

6.154340277778X10 1 
-3.876159660218X10 2 

2.730071323411X10 3 



H = 0.5 
2.000000000000 
3.000000000000 

-6.833333333333 
4.075451388889 xlO 1 

-2.760344845403 xlO 2 
2.011538052494x10 s 



2.000000000000 

3.000000000000 
-9.833333333333 

5.498784722222X10 1 
-3.678279451885X10 2 

2.664691694041x10 s 



2.000000000000 

3.000000000000 
-9.833333333333 

4.683784722222 xlO 1 
-2.870075748181 xlO 2 

2.110416878518x10 s 



3.000000000000 

1 4.666666666667 

2 -9.740740740741 

3 2.117987280210X10 1 

4 -7.103051896801 xlO 1 

5 2.054506982964 xlO 2 

6 -8.384073194275 xlO 2 



3.000000000000 

3.333333333333 
-4.259259259259 

9.052263374486 
-2.746737478771 xlO 1 

9.091443816105X10 1 



ti = l 
3.000000000000 
2.000000000000 

-1.740740740741 
5.541189674523 

-1.941553551960X10 1 
6.590006451270X10 1 



3.000000000000 

2.000000000000 
-3.518518518519 

9.760905349794 
-3.281454899895X10 1 

1.138794166995X10 2 



3.000000000000 

2.000000000000 
-3.518518518519 

6.595922184811 
-2.157173735275X10 1 

7.209915733050 xlO 1 



5.000000000000 

1 2.800000000000 

2 -3.229333333333 

3 3.119580472860 

4 -4.249661491769 

5 9.281773558871 

6 -2.667077094545 xlO 1 



5.000000000000 

2.000000000000 
-1.277333333333 

8.464108513709 x 10" 1 
-1.098495894429 

1.580697986782 



(i = 2 
5.000000000000 
1.200000000000 

-2.213333333333x10- 
4.398750760351 xlO" 

-8.941846899782x10" 
1.236101545345 



5.000000000000 
1.200000000000 
1 -1.181333333333 
1 1.428440057720 
1 -2.414224882476 
4.024451851449 



5.000000000000 

1.200000000000 
-1.181333333333 

9.870369807970X10- 1 
-1.183255224582 

2.130785627667 



2.100000000000X10 1 

1 6.666666666667X10- 1 

2 -5.901090594968X10- 1 

3 -2.582504019689X10- 3 

4 2.457235107452X10" 1 

5 6.900777782581 xlO" 2 

6 -1.752056769026X10- 1 



2.100000000000X10 1 
4.761904761905X10- 1 
-3.564409890941 x 10" 1 
-1.053385459541X10- 2 
1.646791697513X10- 1 
1.213278624644 xlO" 2 



(i = 10 
2.100000000000X10 1 
2.857142857143X10" 1 

-1.625094482237X10" 1 
1.090172643379X10- 2 
7.346442898971 xlO" 2 

-7.752772094654X10- 3 



2.100000000000X10 1 
2.857142857143X10- 1 
-3.620559334845X10- 1 
3.819111465784X10- 2 
1.373283015502X10- 1 
5.188958501563X10- 3 



2.100000000000X10 1 
2.857142857143 xlO" 1 

-3.620559334845X10- 1 
1.199822627856xl0~ 2 
1.747678732566X10- 1 

-8.047408934744 xlO" 3 
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